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Abstract: Two Nal (85% ^^^I and 15% ^'^^I) targets were exposed to a beam of 660-MeV protons. Cross sections for formation 
of 76 residual nuclei were obtained by the induced activity method. The results are compared with other experimental data on 
^^'^l and theoretical calculations by eleven models contained in the codes LAHET3 (using the Bertini+Dresner, ISABEL+Dresner, 
INCL+Dresner, and INCL+ABLA options), CASCADE, CEM95, CEM2K, LAqGSM+GEM2, CEM2k+GEM2, LAqCSM+GEMINI, and CEM2k+GEMINI. 
Most of the models describe spallation products with masses close to the target reasonably well while the reliability of the 
codes differs greatly in the deep spallation and fission/fragmentation regions. The difficulties in describing products with 
A=40-80 by all of the codes tested here except for CEM2k+GEMINI and LAQGSM+GEMINI is related to the neglect of fission (and 
fragmentation) processes for targets as light as ^^^I. 

Introduction 

This work is a continuation of T', where results of determining cross sections for the yields of residual nuclei 
from the ^^^Np and ^''^Am targets exposed to an extracted 660-MeV proton beam of the JINR DNLP Phasotron 
were reported. The first experiments to study transmutation of ^^^I {T1/2 = 1.57 x lO*" y) and some transuranic 
isotopes were carried out in 1996 with heavy-ion beams of the JINR LHE Synchrophasotron using relativistic protons 
of 3.7 GeV [31. Preliminary results of our measured product cross sections from ^^^I targets irradiated with 660 MeV 
protons were published earlier Here, we present final experimental cross sections for our 85% ^^^I and 15% ^^^I 
target irradiated by 660 MeV protons. Using previous measurements on ^^^I targets generally performed at energies 
below 660 MeV analyzed by Molodo and Holzbach HI and reducing them to our proton energy of 660 MeV by linear 
interpolation between energies 600 and 800 MeV, we estimate here experimental cross sections for the target ^^^I at 
660 MeV. We analyze our measurements with eleven models incorporated into several transport codes used recently 
in different applications. Details on our measurement and on models used here to analyze the data may be found in 
Ref. 5, 

Results and Discussion 

Table ^ presents our data: Column 1 lists the isotopes (ground or metastable states) for which we obtained 
reliable production cross sections, columns 2 and 3 show their half-lives (T1/2), decay and cross section types, 
respectively, column 4 presents the measured cross sections for our 85% ^^^I and 15% ^^^I target, column 5 shows 
cross sections for ^^'^l from and column 6 shows our deduced experimental cross sections for ^^^I. 



Table 1: Experimental results 
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Table 1: Experimental results (continued) 
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Table 1: Experimental results (continued) 
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1 1 Sm 
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1 18 J 


13.7 m 


C(/3+,EC) 


3.3(4) 
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16.05 h 


C(EC,/3+) 
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119 TTL ' 


4.69 d 
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154.0 d 


I(IT,EC) 


13.3(16) 


16.9(17) 


12.7(13) 




2.70 d 


CC/3- EC 3+) 




4 3('5'l 


8.3(10) 


123j 


13.2 h 


C(EC) 


25.0(24) 


26.1(30) 


24.6(30) 


124j 


4.18 d 


I(EC,/3+) 


26.3(30) 


32.3(40) 


25.3(35) 




12.4 d 




0.54(5) 






126j 


13.02 d 


I(IT) 


35(5) 


66.0(80) 


29.5(20) 


i27Xe 


36.46 d 


C(EC) 


3.8(4) 


1.42(20) 


4.2(5) 


128j 


25.0 m 


I(/3-,EC) 


31(5) 







We analyzed our measurements with the LAHET3 version IB| of the transport code LAHET |7j using the Bertini 
[H] and ISABEL P| intranuclear-cascade models (INC) merged with the Dresner evaporation model jJOl and the 
Atchison fission model (RAL) and using the code INCL by Cugnon et al. merged in LAHET3 with the 
ABLA and with Dresner JOI (+ Atchison JI]) evaporation (+ fission) models, with the Dubna transport code 
CASCADE |14j. with versions of the Cascade-Exciton Model (CEM) ^ as realized in the codes CEM95 ^ and 
CEM2k [TTj, with CEM2k merged [111-1101 with the Generalized Evaporation/fission Model code GEM2 by Furihata 
,21) . as well as with the Los Alamos version of the Quark-Gluon String Model code LAQGSM [22] merged |18 | -|20 | 
with GEM2 and with versions of the CEM2k and LAQGSM codes both merged with the sequential-binary- 
decay code GEMINI by Charity [23] . Most of these models are widely used to study reaction on heavier nuclei (see, 
e.g., |^|2S| and references therein). 

As we have done previously (see, e.g., [^[2Fj), we choose here two quantitative criteria to judge how well our 
data are described by different models; namely, the ratio of calculated cross section for the production of a given 
isotope to its measured values a'^"-^ /a'^^^ as a function of the mass number of products (Fig. 1), and the mean 
simulated-to-experimental data ratio (Tables 2 and 3) 

(F) = ioV<(i°g[°"°'/'^="''])^> ^ (1) 

with its standard deviation : 

5 ((F)) = ioV((l'°s('^"'/'^''"'')l-'°s«-^»)'') . (2) 

For such a comparison, out of all the 76 cross sections measured in this work, only 48 were selected to satisfy 
some rules based on appreciation of the physical principles realized in the models. For instance, if only an isomer 
or only the ground state of a nuclide with a relatively long-lived isomer was measured, such nuclides were excluded 
from the comparison, but if both were measured separately, their sum was compared with calculations. Such rules 
are essentially similar to those used by Titarenko et al. |24[ |2S] . 

To understand how different models describe nuclides produced in the spallation and fission or fragmentation 
regions, we divided all 48 measured nuclides included in the comparison into two groups, spallation [A > 95) and 
fission/fragmentation [A < 95). Table 2 shows values of {F) and S {{F)) for all compared products (both spallation 
and fission/fragmentation), while Tab. 3 shows such results only for spallation; N is the total number of comparisons, 
-^30% — the number of comparisons in which calculated and measured values differ by not more than 30 %, A''2.o — 
where the difference was not more than a factor of two. 

We note that the codes CEM95 and CEM2k |T2] consider only competition between evaporation and fission 
of excited compound nuclei and calculate the fission cross sections for a nuclear reaction on a heavy nucleus, but 
do not calculate the fission fragments, as they do not contain a fission model. The Bertini [S] and ISABEL [Sj 
INC's are used in our calculations with the default options of LAHET3 for evaporation/fission models; they consider 
evaporation with the Dresner code jJOj and a possible fission of heavy compound nuclei using the Atchison RAL 
fission model [2], but only if they are heavy enough {Z > 71), i.e., do not consider fission for such light targets 
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Figure 1: Comparison of experimental and theoretical cross sections for our 15% ^^^I + 85% ^^^I target. 



as 1291. CEM2k+GEM2 and LAQGSM+GEM2 consider fission using the GEM2 model jH] of only heavy nuclei, 
with Z > 65, i.e., also not considering fission of our target. Similarly, INCL+ABLA ^JE] also does not consider 
fission for L Only the code GEMINI by Charity [53] merged with CEM2k and LAQGSM considers fission (via 
sequential binary decays) of practically all nuclei, and provides fission products from our reactions. This is why 
CEM2k+GEMINI and LAQGSM+GEMINI agree better than all the other models tested here with experimental 
data for this reaction, especially in the A = 40-80 mass region. 

Newer calculations have shown that it is possible to extend the fission model of GEM2 so that it describes 
also fission of light nuclei, like ^^^I, and gives with CEM2k+GEM2 and LAQGSM+GEM2 for our reactions results 
very similar (even a little better) to the ones provided by GEMINI. For this, it is necessary to fit the ratio of the level- 
density parameters for the fission and evaporation channels, a//a„. We think that it is possible to extend in a similar 
way also the Atchison fission model ^I] and the ABLA evaporation/fission model .13_ to describe fission of Iodine 
also with the BertiniH-Dresner/ Atchison, ISABEL+Dresner/Atchison, INCL+Dresner/Atchison, and INCL+ABLA 
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Table 2: Comparison for all 48 (spallation and fission/fragmentation) selected isotopes 



Model 


N/N3o%/N2.o 


(^> 


^((^>) 


Bertini+Dresner 


43/ 6/23 


6.18 


5.28 


ISABEL+Dresner 


40/ 6/18 


5.70 


4.59 


INCL+Dresner 


36/13/21 


4.10 


3.12 


INCL+ABLA 


36/10/19 


9.95 


6.41 


CASCADE 


48/ 9/16 


11.36 


5.20 


CEM95 


46/ 9/22 


5.15 


3.27 


CEM2k 


40/15/29 


6.50 


6.62 


LAQGSM+GEM2 


42/12/23 


14.37 


12.78 


CEM2k+GEM2 


43/11/31 


11.08 


10.80 


LAQGSM+GEMINI 


48/17/32 


4.21 


3.35 


CEM2k+GEMINI 


48/12/29 


2.78 


2.11 



Table 3: Comparison for only 26 selected spallation isotopes with A > 95 



Model 


N/N3o%/N2.o 


(^^> 


sm) 


Bertini+Dresner 


26/ 6/20 


1.86 


1.47 


ISABEL+Dresner 


26/ 6/17 


1.91 


1.48 


INCL+Dresner 


26/13/21 


2.07 


1.89 


INCL+ABLA 


26/10/19 


2.46 


2.07 


CASCADE 


26/ 9/15 


3.79 


3.02 


CEM95 


26/ 8/20 


1.93 


1.51 


CEM2k 


26/12/22 


1.72 


1.51 


LAQGSM+GEM2 


26/12/21 


1.84 


1.61 


CEM2k+GEM2 


26/ 9/22 


1.84 


1.55 


LAQGSM+GEMINI 


26/16/23 


1.55 


1.40 


CEM2k+GEMINI 


26/ 9/20 


1.75 


1.49 



options of LAIIET3; the same is true for the Dubna code CASCADE. Nevertheless, we are not too optimistic about 
the predictive power of such extended versions of these codes as they do not yet contain reliable models for fission 
barriers of light nuclei. 

To make the situation even more intricate, we note that when we merge CEM2k+GEM2 and 

LAQGSM+GEM2 with the Statistical Multifragmentation Model by Botvina et al. [37|, it is possible to describe 
these reactions and get results very similar to the ones predicted by CEM2k+GEMINI and LAQGSM+GEMINI 
without extending the fission model of GEM2, i.e., considering only INC, preequilibrium, evaporation, and multi- 
fragmentation processes, but not fission of ^^^I. We will discuss these results in a future publication. Here, we note 
that it is impossible to make a correct choice between fission and fragmentation reaction mechanisms involved in 
our p+^^^I reaction by comparing theoretical results with only our (or other similar) data; addressing this question 
would require analysis of two- or multi-particle correlation measurements. 

From Fig. 1 and Tabs. 2 and 3 we see that the agreement of different models with our data varies quite a bit. 
We find that most of the codes are fairly reliable in predicting cross sections for nuclides not too far away in mass 
from the target, but differ greatly in the deep spallation, fission, and/or fragmentation regions. We conclude that 
none of the codes tested here is able to reproduce well all of our data and all of them need to be further improved; 
development of a better universal evaporation/fission model should be of a high priority. 
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